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ABSTRACT 


The occurrence of graphite as a common accessory mineral in meteorites 
and in terrestrial metamorphic and igneous rocks gives particular importance 
to the study of equilibrium between graphite and a coexisting gas phase. 


Using a simplified model in which T, P g as > an ^ fg are independently specified 
for the system C-H-0 , values of P , P „ , P u ^ , P u , and P_ in a gas 

OU 2 CaJ H^U ri2 

phase in equilibrium with graphite have been calculated for a wide range of 

geologically possible conditions by use of a high-speed computer. The 

numerical results support the following general conclusions: (1) The 

assumption that P = P + P is significantly in error for many 
gas H 2 O 00 2 

graphite-bearing mineral assemblages. (2) Methane, CH^, may be a significant 
to dominant constituent of the gas phase in many possible geological environ- 
ments involving moderate reduction; in particular, the occurrence of graphite 
with reduced minerals such as fayalite, wllstite, and iron is indicative of 
a methane-rich gas phase. (3) Under metamorphic conditions, pure water is 
not stable with graphite, but graphite can coexist with a gas phase rich 
in CC^. (4) Original graphite in a sediment will stabilize increasingly 
reduced mineral assemblages during progressive thermal metamorphism. (5) 

The presence or absence of even small amounts of graphite can explain P 

U 2 

gradients observed over short distances or between adjacent layers in 
metamorphic rocks. (6) It is possible that the terrestrial atmosphere 
could have evolved by conversion of original methane to water and 00 2 by 
reaction with graphite and other accessory minerals within the primordial earth 
at temperatures of 600° to 1000°G. Material requirements for such conversion 


are not unreasonable, and the process itself is consistent with many proposed 
models for the origin of the earth. 
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Symbol s 


T 

P 


gas 

P. 

1 

f . 


L 



K. 

1 


log N 
In N 


Temperature in degrees Kelvin. 

Total pressure of a gas phase. 

Partial pressure of component i in the gas phase. 

Fugacity of component i in the gas phase. 

Fugacity coefficient of component i in the gas phase, defined by 

Vi * V p i- 

Equilibrium constant of the ith reaction, defined by AG° = -RT In K 
A function of T only when written using fugacities of volatile 
components . 

Logarithm of N to the base 10. 

Logarithm of N to the base e. 


Abbreviations for solid phases 


H 

hematite 

Pe 2°3 

M 

magnetite 

Fe 3°4 

W 

wtistite 

Fe 0 

1 -x 

I 

iron 

Fe 

N 

nickel 

Ni 

NO 

nickel oxide 
(bunsenite) 

NiO 

Q 

quartz 

Si0 2 

F 

f ayal ite 

Fe^iO 

Gr 

graphite 

C 


Abbreviations for solid-phase oxygen buffer assemblages (Eugster and Wones, 


(1962). 





HM 

hematite + magnetite 

NNO 

nickel 

+ bunsenite 


MW 

magnetite + wllstite 

QFM 

quartz 

+ f ayal ite + 

magnetite 

WI 

wtistite + iron 

QFI 

quartz 

+ f ayal ite + 

iron 


MI 


magnetite + iron 


INTRODUCTION 


The system C-H-0 is of extreme importance in consideration of any 
fundamental problems concerning the origin and development of the earth 
and the solar system. Part of the importance is due to the fact that 
the three elements together comprise about 63 per cent by weight of the 
matter in the universe (Allen, 1963, p. 30-33), of which almost all is 
hydrogen. If th£ inert gases helium and argon are excluded, then carbon, 
hydrogen, and oxygen make up about 99 per cent. In addition, the system 
contains such biologically important compounds as water, oxygen, and 
carbon dioxide, which are of prime importance in speculations about the 
origin of the terrestrial oceans and atmosphere and the origin of life. 

Calculations of equilibrium in the system C-H-0 are complicated by 
the large number of organic and inorganic compounds which the system 
contains; high-speed computers have recently been employed to solve the 
numerous equilibria involved. In previous studies, the elemental composition 
of the system has been arbitrarily specified in order to supply sufficient 
constraints for solution of the equations applied to a single ideal fluid 
phase. Suess (1962) has presented calculations involving only the compounds 
H 2 0 , H 2 , C0 2 , 00, and CH^ f calculating partial pressures as functions of 
the composition of the system. A more general calculation involving additional 
elemental components (Dayhoff et al . , 1964) produced additional equilibrium 
concentrations for a number of organic compounds of biological importance. 

Such calculations, in which total composition of the system is initially 
specified, yield considerable information, but the data are not directly 
applicable to a number of important geological problems. In geological 
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processes such as the metamorphism of sediments, it is not generally possible 
to specify the total composition of the system or to assume that the com- 
position has remained constant. Commonly, the composition of the system 
will change during such processes by gain or loss of both volatile and 
involatile components. Thus, the chemical potentials of the components 
become more important than the original compositions in defining the behavior 
of the system at any given time (Thompson, 1955; Korzhinskii, 1959; 

Zen , 1963). 


Further, in geological investigations, it is generally impossible to 

obtain a sample of the original fluid phase. One is thus compelled to study 

the compositions of the solid minerals in order to estimate the nature of 

the fluid phase with which they were equilibrated. 

The present study considers equilibrium between graphite and a gas 

phase in the system C-H-0 over a range of geologically reasonable temperatures 

and gas pressures. It is an extension of work carried out on the system 

C— 0 (French and Eugster, 1965), which considered the effect of graphite on 

P n in geological assemblages. The present study consists of three stages. 
u 2 

First, compositions of the gas phase in equilibrium with graphite are 

calculated by use of a high-speed computer. Second, independent estimation 

of P n for assemblages of iron-bearing minerals is applied to the data to 
u 2 

estimate the composition of the coexisting gas phase. Finally, the data are 
applied to a simplified model in order to suggest a possible method for 


evolution of the terrestrial atmosphere. 
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THEORETICAL CONSTRAINTS ON THE EQUILIBRIUM CALCULATIONS 

The calculations described here have been obtained from a model which 
incorporates the following assumptions: (1) that continuous equilibrium is 

established between crystalline graphite and a gas phase whose composition 
lies in the system C-H-0 ; (2) that the significant species in the gas phase 
are 00 2 , 00, ^0 , H^ , CH^, and 0 2 ; equilibria involving more complex 
organic compounds are not considered; (3) that the gas phase exhibits 
ideal behavior under all conditions. The implications of these assumptions 
are discussed in turn. 

Equil ibrium between graphite and gas 

In the ternary system C-H-0, under suitably reducing conditions, solid 
carbon (graphite at low pressures) can form over a wide range of temperatures 
and gas pressures. The presence of graphite in the assemblage reduces the 
number of variables which must be specified to fix the composition of the 
gas phase. The two-phase assemblage, graphite + gas, may be shown by the 
Gibbs Phase Rule to be trivariant, since F=G-P + 2 = 3 (French and 
Eugster, 1965, p. 1537). 

Graphite, or amorphous carbonaceous material, is a common trace con- 
stituent of numerous sediments and metamorphic rock^ and has occasionally 
been observed in terrestrial igneous rocks and in meteorites. The effects 

of such material on controlling P have been previously suggested (Miyashiro, 

°2 

1964; Mueller, 1964; French, 1964). The assumption that graphite is present 
in natural assemblages is therefote reasonable. The calculations discussed 
are strictly applicable only to crystalline graphite and not to amorphous 
organic material which may originally be present in sediments and low-grade 
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metamorphic rocks. However, such material is observed to change to graphite 
with increasing metamorphism (Quinn and Glass, 1958; French, 1964). 

The model also requires continuous equilibration with a gas phase of 
the proper composition. It does not apply to systems in which a large 
part of the gas phase is composed of additional elements (e.g., N, S, or 
Cl) or to systems which are so dry and impermeable that no gas phase can 
exist. Processes occuring only in the subsolidus, e.g., exsolution of 
graphite from an iron-carbon alloy, are likewise excluded from consideration. 
Consideratio n of inorganic compounds only 

It is assumed here that the compounds OOj, OD , H^O , CH^, and 0^ 

constitute the total pressure of the gas phase, since these are the major 
volatile species involved in reactions between minerals. A detailed calcu- 
lation of equilibrium compositions for a phase containing additional organic 
compounds (Dayhoff et al . , 1964) provides support for this assumption by 
indicating, (1) that the compounds listed do compose almost the entire gas 
pressure, and (2) that the partial pressure of CH^ greatly exceeds that of 
the more complex organic compounds. It is therefore believed that no 
great discrepancy is introduced by this assumption. 

Ideal Behavior of the Gas Phase 

The assumption that the gas phase behaves ideally has been incorporated 
into earlier calculations involving the system C-H-0 (Krauskopf , 1959; 

Suess, 1962; Holland, 1962; Dayhoff et al . , 1964). This assumption probably 
introduces the largest uncertainty into the numerical results obtained in 
the present study. As discussed below, the effect of this assumption is 
difficult to estimate, but it is believed that the present results are 
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satisfactory for seraiquantitat ive considerations and that none of the general 
conclusions is invalidated. 

There is presently little experimental data on the behavior of 
mixtures of the gas species at high pressures, but available data do indicate 
that the behavior will be nonideal at low temperatures. An immiscibil ity 
gap exists in the system R^O-CO^ (Malinin, 1959; Todheide and Franck, 1963; 
Takenouchi and Kennedy, 1964); although there are slight differences between 
the experimental results, it appears that the two-phase region lies at 
temperatures below the critical point of water (373°C) for values of P 

gas 

less than 2000 b. Below 400 b. total gas pressure, the two-phase region 
apparently lies below 350°C. The effect of other compounds on the location 
of the two-phase region is uncertain; large amounts of methane might move 
the critical surface to even low temperatures. 

Preliminary experimental data (H. J. Greenwood, personal communication) 
indicate that the system l^O-CK^ exhibits nearly ideal behavior at higher 
temperatures. Nevertheless, it seems clear that numerical results based 
on an ideal-gas model are only speculative below 30G°C; only the results 
given here for temperatures above 600°K (327°C) are considered as reasonable 
approximations to the actual behavior of the system. 

Although there are no data available for fugacity coefficients of the 
various species in mixtures, values of y for individual species have been 
compiled as functions of temperature and pressure for 00 2 (Kennedy, 1954; 
Robie , 1962), 00 (Newton, 1935), H 2 0 (Anderson, 1964), H 2 (Shaw and Wones, 
1964), and CH^ (Douslin et al . , 1964). A closer approximation to the actual 
behavior of the gas phase could be made by assuming that values of y ^ in 
the mixture are the same as for the pure species at the same temperature and 
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pressure (Lewis and Randall's Rule; see Denbigh, 1957), Even this approxi- 
mation may not be valid at extremely high pressures. 

In view of the uncertainty in for the gas species, it seems unwarranted 
to attempt to compensate for nonideal behavior in the present calculations. 

An estimate of the uncertainty introduced by assuming ideal behavior is 
carried out below by performing parallel derivations for ideal and nonideal 


gas phases. 
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CALCULATIONS OF EQUILIBRIUM BETWEEN GRAPHITE AND GAS 

In a system which conforms to the assumption listed, the various 
possible equilibria between graphite and the gas phase are specified by the 
equilibrium constants of the following four independent equations (French 
and Eugster, 1965, p. 1537): 


(1) C + o 2 = co 2 



_ Y 00 2 P C02 


0. 


, (2) C + 1/2 0 2 * 00 
(3) C + 2 H 2 = CH 4 


„ _ f oo _ Y oo p oo , 
2 ' ' % 2 r 0 2 )% 

_ £ CH ll _ Y CH 4 p ch 4 

3 ' <v s " ^Tv 5 


(4) H 2 + 1/2 0 2 = H 2 0 


K = h 2 ° ,1 

k 4 f„ (f n fi 
H 0 2 


y h 2 o p h 2 o 

y h 2 p h 2 (y 0 2 p o 2 > % 


Combining (1-4) with the condition that Pg as must equal the sum of all 
P^ , we obtain five equations involving six variables, Pqq^’ P qq ’ P H 0’ P H 2 ’ 

P CH,.' and P 0 • 


K 


(la) P 


1 Yo 2 P 0' 


00 , 


f CJ0, 


!i_ f 

^ 00 2 °2 


(2a) P 


CO 


K 2 (y Q 2 P o/ 
Y 00 


K 2 Is 

(f y* 

Y 00 W 


8 


K 3 ( Th P h) 


k 4 Wh p h 2 )( y 0 P 0 / 


V 


p 00 2 + p oo + p h 2 o + p h 2 + P CH 4 + p o 2 


Thus, for any specified T and P , the values of any five P^ are 
explicitly determined if the sixth is independently specified. This result 
agrees with that of the Gibbs Phase Rule; the system graphite-gas is trivariant 
and becomes univariant if any two variables (e.g. , temperature and total 
gas pressure) are specified. 

Theoretically, one can specify any P^ for calculation of the others. 

In practice, in systems involving both natural minerals and synthetic 

phases, it is convenient to specify P (or f ) as the independent 

°2 2 

variable. Numerous mineral assemblages which contain iron- and manganese- 
bearing phases are related by oxidation-reduction reactions which are 
conveniently written in terms of Pq^. Changes in Pq^ will be reflected by 
resultant changes in the mineral assemblage (Eugster, 1959; Eugster and 
Wones, 1962; Buddington and Lindsley, 1964). Conversely, it may be possible 
to deduce approximate values of Pq^ from the mineral assemblage and thereby 
calculate partial pressures of other species. 

The subsequent calculations are therefore carried out by specifying 


T, P , and f q (specifying f Q removes the necessity of considering Yo, 
* 2 * 

in the equations). 
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Equations (la) and (2a) demonstrate that values of P and P 

w 2 W 

depend only on temperature and and are independent of Pg as (except for 
the slight effect of P on y . ) . The values of P and P„ can thus be 

g3.S L UL) 2 CX) 

calculated directly and substituted in (5), as: 


<5a) P CH 4 + P H 2 + P H 2 0 P gas ~ (P 00 2 + P 00 + P 0 2 ) 


The value of Pq^ will be under 10 b. in most geological assemblages, 

particularly if graphite is also present (Eugster, 1959; French and Eugster, 

1965); the term can therefore be neglected. 

The lefthand side of (5a) can be written in terms of P„ and known 

H 2 


quantities by substituting (3a) and (4a). 




+ p H 2 + 


K 4 (Y H, P H > 


(fnJ 


(p oo 2 + P 00 ) 


This expression is rearranged to obtain a quadratic equation in which Pj 
is the only unknown. 


(P u ) 


2 f K 3^H 






Solving, we obtain, 


- K 4VH 2 (f 0 2 ) 


10 



Values of 

PH 2 obtained from the positive root of (7) can then be 
substituted in (3a) and (4a) to derive, respectively, P^^ and P^o ' ‘ ^ l?r an - 
ideal gas, all P^ - f^.and all “-1. Under these conditions, (7) becomes: 


( 8 ) 


-[V £ ‘ 


h 2 


o/ + 1 . 


2K- 


-.2 


{ [yy* * l ] * 4 


S [ P g as - (P °02 * W] f 


2 K„ 


Comparison of (7) and (8) indicates that the difference between values 
of obtained by each method will probably be less than an order of magnitude. 

Available data on the fugacity coefficients of the pure species indicate 
that, above 500 °K, values of y ^ lie within the range 0.7 - 3.0. This variability 
will be reduced because the fugacity coefficients appear in (7) as ratios, 
and values of will tend to be similar for the same temperature and 
pressure (Newton, 1935). 

It is therefore believed that the errors introduced into the calculated 
by postulating an ideal gas phase will be less than an order of magnitude 
and probably closer to a factor of two. Such a discrepancy is not believed 
sufficient to change the general conclusions obtained and may be minimized 
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by presenting results in terms of log P^. Further experimental work 
would be necessary to obtain more accurate values of y^ for refinement of 
the results. One possible refinement would involve using values of y ^ 
for the pure species in accordance with Lewis and Randall’s Rule, but the 
small anticipated correction does not seem to justify introduction of the 
additional complications into the present study. 

Accordingly, calculations of are carried out for an ideal gas, using 

values of T, P g as > an< 3 selected to cover a wide range of possible 

geological environments. Calculations were carried out on an IBM 7090/7094 

system, using a program written in FORTRAN IV language by Mr. Arthur Smith 

of the Theoretical Division, NASA, Goddard Space Flight Center. The 

values of were obtained from JANAF Tables and did not differ significantly 

from those of earlier compilations (Wagman et al., 1945). The calculations 

covered a temperature range of 300° to 1500°K, a pressure range of 1 bar 

to 100 kilobars, and values of fQ^ approximating the interval between the 

quartz-fayal ite-iron and hematite-magnetite buffers (Eugster and Wones, 

1962). The results presented in the text figures were checked by substituting 

the values of P^ into (1-5) to calculate BL 

■^Between 298° and 1200°K, graphite becomes metastable with respect to 

diamond at pressures above about 16 to 40 kb, respectively (Bundy et al . , 

1961). The transformation introduces a slight additional uncertainty into 

the numerical results obtained for values of P above these limits. No 

gas 

correction for this effect has been made, since the probable nonideality 
of the gas phase at such high pressures will introduce a greater error into 
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the results than that caused by the slight difference in AG°, 0.5 to 1.7 
kcal./mole, involved in the diamond-graphite transformation. Furthermore, 
the following discussion is concerned almost exclusively with the region 
within which graphite is stable. 

No correction has been made for the effect of pressure on the 

equilibrium constants of the heterogeneous reactions (1-3), in which the 

solid phase graphite participates. Because of the small molar volume of 

graphite, the change in log for these reactions with pressure is only 

-0.028(P -1 )/T (French and Eugster, 1965, p. 1534). For extreme conditions 

gas 

(Pgag = 10 kb., T = 600°K) , the change in log is less than -0.5; for 
pressures up to 2 kb., the difference will be less than 0.1, and the 
resultant uncertainties do not appear to exceed those introduced by the 
assumption of ideal behavior. 
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NUMERICAL RESULTS AND SOME GEOLOGICAL IMPLICATIONS 


Composition of the Gas Phase as a Function of f^ 

U 2 

The figures and discussion which follow represent only selected samples 
of the large amount of numerical data obtained through use of a high-speed 
computer. The choice is relatively arbitrary, but a selection has been 
made which is intended to apply chiefly to terrestrial geological processes 
under conditions for which the original assumptions remain reasonably valid. 
Figures (1-3) present values of log P. as a function of log f n 

i u 2 

for specified values of T and P . Numerical results for 600 K (327°C) 

gas 

are shown in Figure 1 for total gas pressures of 1 bar, 1000 bars, and 10 

kilobars. This temperature is reasonable for low-grade metamorphism 

(approximately the greenschist facies; Turner and Verhoogen, 1960, p. 534). 

The dotted portions of the curves correspond to values of log f n for which 

2 

the quantity (P^ + P^) becomes greater than P , and for which solutions 
n J 00 2 CD gas’ 

of (8) are not possible. Values of log f^^ for various solid-phase oxygen 
buffer assemblages (Eugster and Wones, 1962) are indicated on the abscissa. 

Similar plots are presented for 900 °K (627°C), a temperature which 
represents probable high-grade metamorphic conditions (Fig. 2) and for 
1200°K (927°C), a reasonable estimate for pyrometamorphism or for some 
igneous phenomena (Fig. 3). 

Several general conclusions appear from the figures. At constant T 
and f Q , higher gas pressures favor a reduced gas phase richer in CH and 
H 2 . Conversely, for a given P ^ and f 02> increasing temperatures produce 
a gas phase which becomes more oxidized and 1^0 and C0 2 appear at the 
expense of CH^. 
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Throughout the temperature range considered, H 9 0 exceeds CCU in the 
gas, even though the gas phase may consist chiefly of other species. 

Finally, for a wide range of geologically reasonable conditions, 
corresponding to the lower part of the stability field of magnetite, methane 
(CH^) is a significant to dominant constituent of the gas phase, particularly 
at lower temperatures (below 600°C) and higher fluid pressures. This 
observation has particular significance for speculations about the fluid 
phase in equilibrium with low- to medium-grade metamorphic rocks and with 
ore minerals in certain hydrothermal deposits. The data presented indicate 
that, at fo 2 values suitable for the appearance of fayalite or wllstite, 
the gas phase is almost entirely methane, particularly below 500°C. The 
amount of methane in the gas phase decreases sharply as T or f increases. 
However, it seems reasonable to state that, where graphite occurs with 
iron-, wllstite-, or f ayal ite-bearing assemblages, the contribution of methane 
to the coexisting gas phase will be significant. At temperatures below 
400°C, the contribution of methane is also important to a gas phase coexisting 
with more oxidized magnetite-graphite assemblages. 

To apply the data to geological environments, it is convenient to 
recast the results presented in Figures (1-3) in order to specify values of 
P^ as a function of temperature. This can be done by determining the value 
of f Q ^ for a given solid-phase oxygen buffer by an equation of the type: 

(9) !og f 02 = - ♦ b + - c(P sas - 12 

(Eugster and Wones , 1962). Thus, for a given buffer and a specified P , 

log f is a function only of T. Values of P. can then be calculated as 
U o 1 

. 2 
unique functions of the temperature. 
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The selection of a specific solid-phase buffer to establish the 
divariant relationship indicated in (9) between T, ^g as » and i m Pl^ es 
the addition of new components to the system C-H-0 and the corresponding 
appearance of additional solid phases in equilibrium with graphite and gas. 

In order that the original trivariant system should become divariant, 

(n + L) solid phases must appear for each additional n components. Thus, 
the addition of Fe produces two-phase buffer assemblages (HM, MI, etc.), 
while addition of both Fe and Si requires three-phase buffers (QFM, QFI). 

In both cases, the new system now has (n + 3) components, (n + 2) solid 
phases, and a gas phase; under these conditions, it is divariant and all 
P^ are explicitly determined by specifying only P^ g and T. 

Because of the rare occurrence of actual buffer assemblages in nature 
(Eugster and Wones , 1962, p. 115-116), it is most useful to use the buffers 
to establish extreme limits of f corresponding to the stability field of 
some index mineral (magnetite, fayalite, wtlstite, etc.); the corresponding 

i 

extreme values of P. can then be calculated. Then, the occurrence of the 

1 i 

mineral with graphite indicates limiting values for composition of the 
coexisting gas phase. 

Although the present discussion has emphasized fg^, t ^ ie technique 

is perfectly general (Eugster and Wones, 1962, p. 95). Additional volatile 

components may be introduced into the present calculations if new f£ can 

be similarly limited by additional solid phases in the system. For example, 

t 

inclusion of fg could be facilitated by introduction of the corresponding 
iron sulfides (Toulmin and Barton, 1964); consideration of f^ might become 
desirable for study of certain nitride-bearing meteorite assemblages (Mason, 
1958; Andersen et al . , 1964; Keil and Anderson, 1965). 
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Most igneous and metamorphic reactions apparently occur at values of 
f 02 within the stability field of magnetite (Eugster and Wones, 1962). 

Figures (4-6) present values of log P. as a function of temperature for three 
buffer combinations within this region: (1) the quartz-f ayal ite-magnetite 

(QEM) buffer (Fig. 4), which fixes the highest fQ^ at which fayalite is 
stable; (2) the combined magnetite-iron (MI) and magnetite-wtlstite (MW) 
buffers, which together specify the lowest f^^ for the existence of magnetite; 
and (3) the combined magnetite-iron (MI) and wtlstlte-iron (WI) buffer's ;■ which 
correspond to the highest fo2 £° r which metallic iron is stable. 

Values of log p^ for the various species in a gas phase coexisting 
with the solid assemblage quartz + fayalite + magnetite + graphite are shown 
in Figure 4 as a function of temperature for P values of 1000 b. and 
10 kb. The effect of a tenfold change in total gas pressure on the relative 
composition of the gas phase is slight; at both pressures, equilibrium is 
not possible above about 550°C. The most important species in the gas phase 
are CH^, CO^ t and 1^0; the latter two become dominant with increasing 
temperature . 

The behavior of the gas phase at P = 1000 b., in equilibrium with 

gas 

the two assemblages, magnetite + iron (below 560°C) and magnetite + wllstite 

(above 560°C) is shown in Figure 5. Methane is dominant over much of the 

temperature range, but the assemblage itself becomes unstable above about 

850°C. The f values for these buffers (Fig. 5) represent minimum values 

for most geological processes. By contrast, at the same values of P , 

gas 

the gas phase coexisting with the more reduced assemblages, magnetite + 
iron (below 560°C) and wllstite + iron (above 560°C) is considerably different 
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(Fig. 6). In this case, equilibrium is possible up to the limit of calculation, 

above 1200°C. The contribution of methane is enhanced, and, at temperatures 

above 450°C, P^ > P^O an ^ Pqq > Pqq • These curves represent the maximum 

f n values for the occurrence of metallic iron, as in meteorites and in 
z 

rare terrestrial rocks. Lower values of f within the iron stability field 
will increase the reduced character of the gas phase. 

Significance of Methane in Equilibrium with Graphi te 

Much speculation about the nature of hydrothermal fluids has considered 
a fluid phase as a mixture of 1^0 and 00^- The present results indicate 

that, in addition, methane may be a significant species under a range of 

geologically possible conditions involving low to moderate temperatures, 
high gas pressures, and moderately low values conditioned by the presence 
of graphite. 

The importance of methane under such conditions is indicated by the 

data presented in Figures (7-8) , which show the percentage of methane in a 

gas phase in equilibrium with graphite and various buffer assemblages. 

Since values of P^ rather than log P^ are plotted, the assumption of ideal 

behavior will cause a proportionately greater uncertainty, but the relative 

relationships should not be affected. At P = 1 b. (Fig. 7), more reducing 

gas ^ 

values of f n favor increased contribution of methane at higher temperatures. 

2 

However, even at values of f n corresponding to the upper part of the 

u 2 

magnetite stability field (approximately the NNO buffer), methane may still 
contribute 30 to 50 per cent of the value of P 

gas 
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Higher gas pressures favor increased contribution of methane at higher 

temperatures; at p g = 1000 b. (Fig. 8 ), methane accounts for more than 

half of the value of P within the stability fields of fayalite and 

gas 

wUstite. Increasing Pg as above 1000 b. will correspondingly increase the 
relative contribution of methane. 

With increasing temperature, methane gradually disappears as ^0 and 00 2 

became dominant. The relative amounts of the latter two species depend on 

P , as indicated in Figure 9, where values for all P. are plotted as a 
gas’ i 

function of T for the QFM and MI-WI buffers. At p = 1 b. , CH /t and H 0 are 

x gas 4 2 

dominant at low temperatures; with increasing temperature, P passes through 

a maximum and C0 o becomes the dominant species. By contrast, at P = 1000 b., 
2 gas 

(Fig. 10), both H 2 and 00 are negligible. With increasing temperature, P^ Q 
and P^ increase continually as P^ decreases, with P u n always greater 


CH 


H 2 0 


than P . Pressures above 1000 b. will increase the preferential production 

o ° 2 

of water with decreasing P , although somewhat higher temperatures will 

CH4 

be required for complete conversion. 

A significant feature of these data is the rapid decrease in P^ with 

OH 4 

increasing temperature. For example, at P g as = 1 b. (Fig. 9), drops from 

0.7 b. to virtually zero between about 225° and 325°C along the QEM buffer; 
even at lower values of f n along the MI-WI buffer, P is reduced from 0.6 b. 

u 2 UH 4 

to zero between 450° and 650°C. Similar but less sudden decreases occur at 

o 

P = 1000 b. (Fig. 10); on both buffers, a temperature interval of only 150 C 
ga s 

is sufficient to reduce P from about 600 b. to zero. Such rapid changes 

CH 4 

in the composition of the gas phase with changing temperature suggest that the 

principal difficulty in estimating gas compositions from natural mineral 

assemblages will be uncertainties in temperature; the effects of uncertainties 

in P will be minor in comparison. The application of these data to geological 
ga s 
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assemblages will thus require very accurate temperature data, since a change of 
only 200° C is sufficient to convert a gas phase composed almost entirely of 
methane into one in which 1^0 and CD^ together constitute more than 80 per 
cent of ^g as * The implications of these results for a model of production of 
the terrestrial atmosphere by reactions within the earth is discussed later. 

The detection and identification of methane in liquid inclusions and in 
some graphite-bearing rocks provides experimental support for the conclusion 
that methane may be a significant constituent of natural hydrothermal fluids. 
Methane has been reported as a minor constituent of liquid inclusions in 
quartz-bearing pegmatites (Elinson and Polykovskii, 1963a, 1963b) and in large 
amounts in inclusions in alkaline igneous rocks, associated with small amounts 
of organic matter (Peter sil T ye et al . , 1961)* At least 5000 ppm of methane 
and minor amounts of other hydrocarbons were evolved from graphite from a 
Precambrian quartz -graphite schist (Hoering and Hart, 1964). Further studies 
of this type would be very desirable. 

Maximum £g ag - Temperature Curves for Graphite - Bearing Assemblages 

In the calculations presented, it was often observed that, under conditions 

of relatively high T and f n , no solutions were possible to the system of 

u 2 

Equations (1-5, 8). This incompatibility arises from the fact that values 

of P m and P are (for the ideal -gas model) dependent only on T and fQ 
2 CD 2 

and completely independent of Pg as * it i- s therefore possible, as T or f^^ 
increases, to generate a condition where the quantity, (Pqq + P qq) exceeds 
the specified value of P g . Such a condition is impossible and the equations 
become insoluble. 


Under conditions of specified land P g as > this limit is reached by 

increasing f Q (Figs. 1-3). In later calculations the use of a buffer to 
2 
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control f Q as a function of T introduces an additional equation involving 

f . Thus, two simultaneous equations must be satisfied for equilibrium 
°2 

to be maintained. For graphite, 


(1) G + 0 2 = 00 2 


■ V £ 0 2 


The solid-phase buffer has an equation of the form 


(10) X + Y = Z + 0, 


K io ■ £ o. 


where X, Y, and Z are solid phases. 

For an ideal gas, P = K f . Further, the system requires that 

w ^ w 2 


P S P_ = K, K = K f. 
gas 00 2 1 10 1 °2 


No solutions will exist for values of T, P and f for which 

°2 


gas 0, 


log P gas * l0g P 00 = l0g K 1 + l0g P - 


0, 


or, 


log p a log p gas - log Kj 


( 11 ) 
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With increasing T or f Q , such a limit occurs below 1500°K for all but 

2 

the more reducing buffers and corresponds to the appearance of disequilibrium 
between graphite and the solid-phase buffer assemblage. In general, the more 
oxidizing the buffer, the lower the temperature at which it becomes unstable 
with respect to graphite. Actual disequilibrium temperatures will be lower 
than suggested by (11) because of the contribution of other species to P 

gas 

Calculation of the exact disequilibrium temperature for different values 

of Pg ag establishes a maximum temperature for equilibrium between graphite and 

each buffer assemblage. These points form a series of univariant curves in 

P -T - f 0 space and may be presented as projections onto the P -T 
gas u 2 

plane (Fig. 11). 

Each curve fixes only the maximum temperature for existence of the buffer 

assemblage with graphite; the phases will be in equilibrium at all temperatures 

below the curve. For instance, the assemblage, hematite + magnetite + 

graphite is stable at all temperatures and pressures to the left of Curve (1) 

(Fig. 11); the value of f in this region is specified by the HM buffer. 

As Curve (1) is reached with increasing temperature, the three phases become 

unstable and one must disappear. Three paths are possible: (1) graphite is 

removed by oxidation, and the assemblage hematite + magnetite continues to 

buffer f at temperatures to the right of the curve; (2) hematite is reduced 
u 2 

to magnetite, in which case magnetite and graphite coexist to the right of the 

curve and do not constitute an fo 2 buffer; (3) magnetite is oxidized to 

hematite, producing the assemblage hematite + graphite, which also does not 

buffer f Q . If graphite remains in the assemblage at temperatures to the right 
2 

of the curve, the buffering ability of the system is destroyed and two 
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assemblages are possible: hematite + graphite at high f 0 and magnetite + 

2 

graphite at low fc^- 

Similar arguments apply to the other curves, particularly to Curve (3), 

which indicates the maximum stability temperature for the assemblage, quartz 

+ fayalite + magnetite + graphite. The four-phase assemblage, stable to the 

left of Curve (3), buffers f n rigorously. To the right of the curve, four 

U 2 

three-phase assemblages are possible: (a) quartz + fayalite + magnetite; 

(b) quartz + fayalite + graphite; (c) quartz + magnetite + graphite; (d) 
fayalite + magnetite + graphite. Only assemblage (a) is a rigorous oxygen 
buffer. 

In general, for each curve, the region to the left denotes equilibrium 

between the solid-phase oxygen buffer, graphite, and the gas phase. As the 

curve is crossed with increasing temperature, the buffer assemblage and 

graphite become unstable, and the character of the system at higher temperatures 

depends on whether or not graphite disappears. If graphite is oxidized, the 

solid-phase buffer still controls f , but another degree of compositional 

U 2 

freedom is added to the gas phase by the disappearance of graphite. If graphite 

remains, the solid phases no longer constitute an oxygen buffer and the value 

of fQ can be specified only within wide limits. 

The curves have considerable application to geological problems, for they 

provide upper temperature limits for the existence of certain mineral assemblages 

in equilibrium with graphite and gas. They have the added advantage of being 

extremely insensitive to changes in P above a moderate pressure. The difficulty 

gas 

in applying the curves to natural assemblages lies in the necessity of 
identifying several solid phases in the same specimen and of demonstrating 
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that they are in fact in equilibrium (see Zen, 1963). In addition, even if 
it can be shown that a given curve has in fact been exceeded, careful 
study is necessary to distinguish between the different petrologic inter- 
pretations deriving from the presence or absence of graphite in the assemblage. 

The fact that most present solid-buffer data involve iron-bearing mine- 
rals has restricted their application to metamorphosed iron formations 
containing some of the required assemblages (Gundersen and Schwartz, 1962; 

Zen, 1963; French, 1965). However, numerous possibilities exist for study 
of graphite and opaque accessory minerals in igneous and metamorphic rocks 
(see, e.g., Chinner, 1960; Buddington andLindsley, 1964). 

Progressive Metamorphism of Graphite - Bearing Sediments 

There has been much consideration of the role of water and 00^ in 
metamorphism, since many of the observed mineral transformations can be 
expressed by reactions involving these volatile components. Many speculations 
about the nature of hydrothermal fluids have presupposed a gas phase in which 
these two species are dominant. Even though other components such as methane 
may assume a significant role under certain conditions, it is worthwhile to 
consider the relationship between q and Pqq^ ^- n a § as ph ase i n equilibrium 
with graphite under conditions corresponding to progressive metamorphism. 

The presence of graphite or amorphous carbon in metamorphosed sediments 

has been proposed to explain two observed characteristics of metamorphosed 

rocks: (1) apparent relative reduction of mineral assemblages with increasing 

metamorphic grade, corresponding to an increase in the Fe^ + /Fe^ + ratio 

(Shaw, 1956; Zen, 1963; Miyashiro , 1964), and (2) preservation of apparently 

steep gradients in P n over short distances, particularly between adjacent 

u 2 

layers of differing composition (Chinner, 1960; Kranck, 1961; Zen, 1963). 
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Application of the present model requires only the development of 
graphite during metamorphism, possibly through recrystallization of original 
organic matter (French, 1964; Hoering and Hart, 1964). Because of the low 
gram-molecular weight of graphite, a small amount will be ample for expensive 
equilibration; such amounts are common in metamorphosed sediments (see, 
e.g., Miyashiro, 1964). 

Construction of q-Pqq plots from the numerical results requires 

only the selection of P m as an independent variable instead of f n . Then, 

2 u 2 

for a given temperature and P , P ' can be plotted as a function of P~. 

gas H 20 r OO 2 

to form a series of isothermal curves (Figs. 12-13). Each point on a curve 
corresponds to a condition in which all P^ (including Pq^) are explicitly 
specified . 

At P = 100 bars (Fig. 12), the P u n -P^ isotherms between 300° 
gas ul) 2 

and 1500°K are a series of asymmetric curves which pass through maxima in 

the region of 5 to 30 b. P^ . Additional curves corresponding to data 

derived from the f 0 -T relations for the solid-phase oxygen buffers have 

2 

been superimposed on the isotherms. 

The oblique line at the top represents the locus of points for the 

condition, P Q + p = 100 b., and the vertical distance between this line 
2 2 

and any point on an isotherm represents the total pressure of other species 

in the gas phase. At low values of P__ , this difference is large for all 

(JO 2 

temperatures and reflects the dominance of methane under conditions of low 

V At higher temperatures, the difference is large for all values of 

P 00 2 because of the increasing contribution of H^ and 00 to the gas phase. 

An increase in P^ from 100 b. to 1000 b. (Fig. 13) produces a slight 

shift in the isotherms toward higher values of P and P^ ; however, the 

H^O 00 ^ 
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general relations are unchanged. Increases in P above 1000 b. will 

gas 

produce relatively smaller shifts. 

These relations have several implications for study of the progressive 

metamorphism of graphite-bearing sediments. First, graphite cannot coexist 

with pure water under the conditions indicated; even at the lowest temperatures, 

the highest water content permitted is about 0.85 P . By contrast, graphite 

gas 

can coexist over a wide temperature range with a water-free atmosphere con- 
taining chiefly CD^. 

Second, for a fixed temperature and P , the values of P„ ^ and P m 

K gas’ H 2 0 00 2 

cannot vary independently but must correspond to points along an isotherm. 

As P increases, P increases at first, passes through a maximum in 

uu 2 h 2 u 

the region where P > P , then decreases to zero. 

00 2 

For the complete temperature range under discussion, low values of 

P (less than 0.05 P ) reflect a methane-rich atmosphere existing under 
oo 2 gas 

conditions of relatively low Pq . With increasing temperature, however, 

the gas phase becomes progressively depleted in H 2 0 and OO^and H 2 and CD 

become important. For this reason, the assumption that P = P^ + P TT ^ 
y gas 00 2 H 2 0 

is significantly in error over a wide range of metamorphic temperatures. 

At 600°K (327°C), the contribution of other components will exceed 10 per 

cent of P for values of P^ ^ 0.4 P . At temperatures above 1000°K 
gas 00 2 . gas y 

o 

(727 C) , a similar error is present for all values of P . 

uj 2 

Comparing the solid-phase buffer curves with the isothermal P -P 

H 2 U 00 2 

curves (Figs. 12-13) indicates that reaction with graphite during rising 
temperature will produce progressively reduced mineral assemblages. For 


instance, at P g = 1000 b. (Fig. 13), temperatures between 100° and 200°C 
are sufficient to convert original hematite to magnetite if graphite is 
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present. Such a mechanism has been suggested to explain the development 
of magnetite from original hematite in slightly metamorphosed sedimentary 
iron formations (Zen, 1963; LaBerge, 1964). More extreme reduction may be 
attained at higher temperatures; the development of fayalite in contact- 
metamorphosed iron formation may be conditioned by the presence of graphite 
(Gundersen and Schwartz, 1962; French, 1965). The rare occurrences of 
wtlstite and native iron in terrestrial rocks (Walenta, 1960; Mel son and 
Switzer, 1965) also appear associated with graphite which, at igneous 
temperatures, could produce extremely reducing conditions. It should be 
added that, when such reduced minerals are observed, the contribution to 
P of species other than H,0 and 00, may be significant. 

A specific model of particular interest is the progressive metamorphism 
of a sediment containing graphite and original interstitial water under 
conditions where the system is either completely closed or open only to 
water (Yoder, 1955; Thompson, 1955). It has been pointed out by earlier 
investigators (Eugster, 1959; Miyashiro, 1964) that, under these initial 
conditions, the resultant gas phase must undergo reduction in order to be 
in equilibrium with the observed mineral assemblages. 

In applying the present calculations to such a system, we may express 

the additional constraint by saying that, whatever the composition of the 

gas phase during metamorphism, the total H/O ratio is fixed at 2/1 for both 

the closed and water-permeable systems. Such arbitrary selection of elemental 

ratios to reduce variance has been applied previously in analogous calculations 

(Suess, 1962; Dayhoff et al . , 1964). This additional restriction makes it 

possible to calculate all P. as explicit functions of T and P only. In 

l gas 

addition to the equilibrium relations (Eqs. 1-5), the requirements of 
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materials balance may be written (Suess, 1962): 


H = 4 P„ + 2 P + 2 P 

ch 4 h 2 h 2 o 


0 = P +P +2P + 2 P 

r H 2 0 r 00 CD 2 0 2 


Since H/0 = 2/1, 


4 P <V 2 P H 2 * 2 p h 2 0 - 2 <p h 2 0 * p oo * 2 p 00 2 * 2 V 


<12) 2 p 00 2 * p ao - 2 P CH 4 - 2 P H 2 * 2 p o 2 = 0 


Equations (1-5) and (12) together constitute a system of six equations 

for the six unknown P^; explicit solutions are therefore possible once T and 

P are fixed. In practice, however, the calculations became involved 
gas F 

because of the appearance of a quartic equation in P . Therefore, a simpler 

H 2 

procedure was adopted by specifying T and f instead. P^^ was calculated 


gas 


by summation at the end of the process and values of P^ were obtained by 
interpolation for any selected P 

Under these conditions, graphite becomes a true oxygen buffer, since 
the value of f 0 ^ maintained in the gas phase depends only on T and P g as * 
Values of f q fixed by graphite in this system for values of Pg as of 1 b. 
and 2000 b. (Fig. 14) are very close to analogous values calculated earlier 
for a gas phase composed only of 00^ and CD (French and Eugster, 1965). 

Both sets of f^ -T curves have shallower slopes than do the curves for the 
iron-bearing solid phase buffers. Thus, in both cases, more reduced mineral 
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assemblages will appear with increasing temperature, in agreement with the 

general conclusions discussed above (Figs. 12-13). 

The interpolated values of P^ are shown in Figure 15 as a function of 

temperature for P = 1000 b. Below about 1000°G, the gas phase consists 
gas 

chiefly of H^O , 00 2 > and CH^; the water content decreases sharply with rising 

temperature. P and P are approximately equal and pass through maxima 

(JJ 2 

where each constitutes about 0.25 P . Above about 1000°G, and GO 

gas 2 

become dominant, although the gas phase still contains significant amounts 

of all five species. The f -T curve for this system intersects progressively 

°2 

reduced buffer curves at slight angles with rising temperature. In this 
idealized case, fayalite would form in rocks of the proper composition at 
about 540°C, wllstite at 850°C, and metallic iron at 1160°G. 

Mineral Assemblages in Meteorites 

Numerous recent studies of the mineralogy and petrology of meteorites 

have attempted to apply mineralogical data to estimate conditions of formation 

of the various classes (see, e.g., Ringwood, 1961; Wood, 1962, 1963; Mueller, 

1964). Considerable speculation has involved estimates of P and possible 

0 2 

theories for evolution of the various types under differing degrees of 
reduction . 

The general relationships between graphite and the composition of the 
coexisting gas phase suggest that the occurrence of graphite as an accessory 
mineral in meteorites is of great significance in attempting to evaluate 
conditions of formation. Graphite has been observed in several classes of 
meteorites, most prominently in the irons, but also in smaller amounts in 
stony meteorites, notably chondrites (Mason, 1958; Krinov, 1960; Ramdohr , 

1962; Andersen ,et al . r 1964; Vdovykin, 1964). 
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In particular, the assemblage graphite + nickel -iron is common and 

widespread in both iron and stony meteorites. The results of the present 

study suggest that this assemblage formed under extremely reducing conditions 

in which methane may have been a dominant component of the coexisting gas 

phase, even at relatively high temperatures. For instance, at P = 1000 b. 

gas 

(Figs. 5, 8), the upper limit of stability of iron + graphite is fixed by 
the MI and WI buffers; at all temperatures below 1000°C, P exceeds 500 b. 

ch 4 

Lower values of within the iron stability field will increase the con- 
°2 

tribution of methane, even at higher temperatures. The reported occurrence 
of quartz and tridymite in several iron-bearing meteorites (see Anders, 1964) 
suggests that at least some meteorites have developed at very low values of 
below the QFI buffer where silica and iron can coexist without formation 

°2 

of fayalite. At these values of P (Fig. 8) the same P could be main- 

u 2 GH 4 

tained at temperatures in excess of 1200°C. 

The use of buffer curves for pure iron in this discussion will not change 

the general conclusions. Even though the nickel -iron in meteorites will be 

in equilibrium with oxide at higher f^ values than pure iron, calculations 

suggest that the effect will not be large for small nickel contents (Olsen, 

1959). Furthermore, the metal phase in at least one class of meteorites, 

the basaltic achondrites, is strongly depleted in nickel (Duke, 1965). 
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A POSSIBLE MECHANISM FOR EVOLUTION OF THE TERRESTRIAL ATMOSPHERE 

The origin and evolution of the atmosphere and oceans are of profound 
importance in consideration of such problems as the origin of the earth and 
solar system, the nature of ancient geological processes, and the origin 
of life on earth. The purpose of the present discussion is not to review 
or evaluate the numerous theories and voluminous literature on the subject 
(for summaries, see Rutten, 1962; Brancazio and Cameron, 1964; Donn et al . , 

1965), but to discuss certain possible processes related to the occurrence 
of graphite and other minerals in the earth. 

Three general methods have been used by earlier investigators to 
estimate the character of the primordial or Precambrian atmosphere existing 
at a time before deposition of most of the presently preserved sedimentary 
record : 

(1) Calculation of the amounts of volatile compounds now present in 
the atmosphere and oceans which must have been produced by processes other 
than rock weathering (Rubey, 1951); such "excess volatiles" must either have 
been present in an original atmosphere or must have been added by gradual 
outgassing from the interior of the earth during geologic time. 

(2) Use of thermodynamic equilibria to determine relative amounts of 
different volatile species in the atmosphere (Rubey, 1951; Holland, 1962, 

1964). Such results are of limited value because they may not consider all 
possible equilibria, and because the equilibrium constants involved are 
strongly temperature-dependent. Temperatures existing in ancient atmospheres 
and within the primordial earth are speculative and the arbitrary selection 

of high temperatures is not justified. Furthermore, if low surface temperatures 
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are postulated, one cannot assume that equilibrium is established, either 
within the atmosphere or between the atmosphere and surface rocks, any more 
than it is at present, 

(3) Study of ancient sedimentary rocks to estimate the character of 
the atmosphere in which they formed. Many studies have cited the occurrence 
of reduced minerals in sediments as evidence that the Precambrian atmosphere 
was reduced with respect to the present one (Rankama , 1958; Holland, 1962; 

Lepp and Goldich, 1964). These studies, while valuable, are hampered by 
two conditions. First, there is no evidence that detrital minerals are 
necessarily equilibrated with the atmosphere during erosion, transport, and 
deposition. Such equilibration is generally not attained at present; magnetite, 
iron-titanium oxides, and ferrous silicates are common detrital accessory 
minerals, although they should oxidize to hematite or iron hydroxides in the 
present atmosphere. Second, the effects of burial, diagenesis, and low- 
grade metamorphism on the mineralogy are difficult to recognize. The sediment 
may be effectively insulated from the atmosphere in a local reducing environ- 
ment soon after deposition. Further, the presence of organic matter or 
graphite may produce reduced mineral assemblages with very slight increases 
in temperature, destroying traces of original equilibrium with the atmosphere. 

It has been generally assumed that a reducing atmosphere lacking free 
oxygen existed at some point in the history of the primordial earth. Such 
an atmosphere would be required for the development and preservation of 
organic compounds and their eventual combination into primitive forms of 
life (see Rubey , 1951; Urey, 1952; Oparin, 1957; Miller and Urey, 1959; 

Holland, 1962; Rutten, 1962). The exact character of this atmosphere is not 
agreed upon; it has variously been proposed that CH^, QO^, and H^O were 


dominant species. 
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Holland (1962) proposed a three-stage model for development of the atmosphere 
through outgassing of the original crust: (1) an atmosphere formed in equil- 

ibrium with metallic iron at 1200°C, composed of CH^ , and minor H^O and 
liquid H^O would also exist on the surface. (2) A later stage containing 
dominant N 2 and minor H 2 0 and 00 2 , produced by equilibrium with magnetite in 
a more oxidized crust; (3) A final stage developed from photo synthetic production 
of oxygen and corresponding to the present atmosphere. 

This model has two serious objections. First, the temperature of 1200°C, 
while not unreasonable, is arbitrary; because of the strong dependence of 
equilibria on temperature, arbitrary selection of temperature strongly 
conditions the gas compositions obtained. Second, the equilibria in the atmosphere 
(at an assumed temperature of 298°K) may not be established, and the calculated 
compositions may not be attained. 

The present discussion is somewhat similar to Holland’s model, and considers 
the problem of converting an original methane-rich atmosphere into water and 
00 2 through equilibria between a gas phase and accessory minerals existing to 
moderate depths within the primordial earth. In the proposed mechanism, 
original methane reacts with graphite and ferric-iron-bearing phases at 
temperatures in excess of 600°G to depths of a few hundred kilometers, producing 
water and C0 2 which are removed to a cool surface to form oceans and atmosphere. 
Such a mechanism does not require low-temperature equilibrium between atmosphere 
and the surface rocks or between the different species within the atmosphere 
itself; the surface temperature can be low enough to prevent re-equilibration 
of water and C0 2 to produce methane again. 

Such a mechanism is independent of any detailed assumptions about the 
composition, thermal history, and structure of the primordial earth. Its 
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operation requires only a source of sufficient original methane and accessory 
minerals, sufficient heating to attain temperatures of 600° to 1000°C at 
moderate depths, and some mechanism for removing the product gases to the 
surface through fissure systems or by volcanism. In the present discussion, 
it is assumed that all of the "excess volatile" water and 00 ^ (Rubey, 1951) 
originated by this mechanism; it will be found that the material requirements 
for even this extreme case are not excessive. Despite the obvious oversimpl if i- 
cations, this mechanism is consistent with a variety of earth models, even 
though it neglects such significant problems as the source of the original 
methane, the time required for the transformation, and the presence or absence 
of a differentiated core in the primordial earth during the process. 

Data from the present study indicate that methane originally in equilibrium 

with graphite will be converted to water and 00 ^ with increasing temperatures; 

at higher values of P g ag » water will form in excess of 00^ (Figs. 9-10). 

At P ^ = 1000 b., temperatures above 600°C are sufficient for complete 

conversion at values of Pq within the stability field of magnetite (QFM 

buffer) (Fig. 10). At minimum values of P within the magnetite field, the 

u 2 

maximum temperature is about 950 C. Increasing P g to 100 kb. increases the 
maximum temperature to about 1250°C. In any case, significantly lower 
temperatures will suffice if water and 00^ are continually removed, allowing the 
reaction to proceed to completion. 

MacDonald (1959) has calculated temperature distributions resulting from 
radioactive heating in a number of primordial earth models as a function of 
time after formation. The numerical results may be subject to later modification 
(Wasserburg et al . , 1964) but it appears that, for numerous possible models, 
temperatures of 800° to 1500°C could be attained at 100 to 300 km. depth within 
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10 years after formation. In fact . MacDonald has suggested (1959, p. 1987- 
1989) that initial temperatures well above these values could have been 
reasonably attained during formation of the earth by accretion. If this were 
the case, methane would have been converted to water and CC> 2 almost simultaneously 
with formation of the earth and would not have existed within*the body for any 
significant length of time. In any case, it appears that the temperatures 
required could have been generated, even though an exact specification of the 
temperature distribution in the body is not possible. 

A wide variety of heterogeneous reactions may occur within the earth 
to convert methane into water and C0 2 through reaction with solid phases; 
the following are some simple examples: 

(13) CH 4 + 2 Fe 3 0 4 + 3 SiC> 2 = 3 Fe 2 SiC> 4 + C + 2 

(14) CH 4 + 2 Fe 3 0 4 + 6 SiC> 2 = 6 FeSi0 3 + C + 2 H 2 0 

(15) CH 4 + 2 Fe 3 0 4 = 6 FeO + C + 2 H O 

(16) 2 CH 4 + Fe 3 0 4 = 3 Fe + 2 C + 4 H O 

Actual reactions will probably involve more complex ferrous-ferric 
silicates such as garnets, amphiboles, and micas, but the stoichiometry will 
be similar to these examples. The subsequent calculations are made using 
magnetite as a reactant and assuming a conversion equation of the form, 

1 CH 4 + 2 Fe 3 0 4 = 2 H 2 0 + C + 3R 


(17) 
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where R designates a reduced phase such as fayalite, ferrosilite, wtistite, 
or iron. The necessary CC> 2 may be formed by reaction of a small amount of 
product water with graphite, i.e., 

(18) C + 2 H 2 0 = C0 2 + 2 H 2 

The volatile products are transported to the surface, where water 
condenses and CO 2 and H 2 are incorporated into the original atmosphere, 
together with an undetermined amount of N 2> If the rate of introduction of 
H 2 is slow, it may escape from the atmosphere without building up a high 
partial pressure (Holland, 1962). 

The amounts of "excess volatile" water and 00 2 which should be produced 
are (Rubey , 1951): 

H 2 0 : 1.66 x 10 24 gm. = 9 . 21 x 10 22 moles 

<X> 2 : 9.1 x 10 22 gm. = 2.07 x 10 21 moles 

From (17) and (18), the materials balance becomes (in moles): 


Reactants Products 


ch 4 

4.81 

x 10 22 

H 2° 

9.21 

X 

io 22 

Fe 3°4 

9.63 

x 10 22 

oo 2 

2.07 

X 

io 21 




H 2 

4.14 

X 

i-H 

CM 
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The amount of methane required, 4.81 x 10 ' moles, would occupy an STP 
12 3 

volume of 1.08 x 10 km. , approximately that of the present earth. If 
this amount had originally existed as an atmosphere, it would have exerted a 
surface pressure of about 140 b. and would have been about 300 km. deep 
(scale height 15 km.). These values are probably excessive and suggest that 
any original methane was probably produced gradually or existed originally 
in the primordial earth. 

Some simple calculations will show that this required methane could have 

been easily contained within a primordial earth as either free gas or organic 

compounds which could yield methane on heating. For discussion, we can postulate 

a primordial earth accreted out of material similar to the carbonaceous 

chondrites, although the chemical composition is not critical to the discussion. 

3 

A body of the same mass as the present earth, with a density of 2.5 gm./cm. , 

would have a radius, R = 8300 km. and a volume, V = 2.39 x 10 1 2 km.^. The 

o o 

23 

required methane, 7.72 x 10 gm. , would be less than 0.00015 of the mass of 
this body. For a primordial earth with a mean porosity of 1 per cent in its 
upper part, a surface pressure of 1 b., a uniform pressure gradient, and a 
constant original temperature of 298°K, the volume necessary to contain the 
required methane can be derived. 

At a depth of h km. from the surface, the lithostatic pressure in bars 
is equal to pgh = 250 h, and the average pressure from the surface to a depth 
h, is P = 125 h. If Boyle’s Law applies to the gas, then V^, the volume of 
rock required to contain the gas, can be calculated from the relationship: 

V STP = °* 01 ^ V h = 1,08 x 1C)12 ktn * 3 
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QTp 11 _ 1 

V L = 0.8 = 8.64 X 10 h 

h a 

For a shell from the surface to a depth h, 

v = v _ 4/3TT (R - h) 3 
ho o 

Solution of the equations indicates that a shell about 35 km. thick would 

contain all required methane. The gas could be contained at a larger depth 

in a correspondingly smaller volume. The exact results are uncertain because 

of the numerous assumptions, but it appears that a primordial earth of any 

reasonable composition might be expected to contain sufficient methane for 

production of the "excess volatiles” by the mechanism suggested. 

It is not unreasonable to assume that the material in the upper region 

of the primordial earth contained at least 1 per cent by weight of magnetite; 

such a condition is met by many terrestrial rocks of widely varying chemical 

composition and also by the carbonaceous chondrites (Mason, 1958). For such 

3 13 

material, each km. of rock contains 5.18 x 10 moles of magnetite and 
11 3 

4.30 x 10 km. of rock will supply the magnetite needed for complete reaction. 
This volume of rock corresponds to a sphere of radius 4680 km., or to a shell 
from the surface of the present earth to a depth of 990 km. (the corresponding 
thickness for the primordial earth postulated above would be about 550 km.). 

It is not necessary to suggest that reaction actually occurred at such 
great depths if some method of convection operated to bring fresh material 
to shallower levels. For example, if reaction can occur to a depth of only 
1 km. , a total of 850 overturn cycles would be needed to supply the required 
magnetite. The number of cycles required decreases approximately linearly as 
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the thickness of the reacting layer is increased; thus, a 20-km. layer requires 

about 40 renewals, and a 100-km. shell only 9. If such convection takes place 

at rates estimated for the present earth (about 1 cm./yr.; Orowan, 1964), 

8 

material could move to the surface from a depth of 1000 km. in about 10 
years. The operation of such a mechanism would not change the total amount 
ot magnetite required, but would allow reaction to be restricted to a narrow 
zone closer to the surface. 

21 

The amount of hydrogen produced by complete reaction, 4.14 x 10 moles 
21 

or 8.34 x 10 grams, would exert a maximum surface pressure of about 1.5 b. 
However, loss of hydrogen into space would tend to keep P at a lower value 

h 2 

(Holland, 1962). 

21 

The total C0 2 produced, 2.07 x 10 moles, would have exerted a surface 
pressure of about 17 b. Precipitation of carbonate sediments would probably 
have occurred long before such a high value was reached (Rubey, 1951) and it 
is questionable whether life could have originated and survived under such a 
high 00 2 pressure. These results suggest a gradual evolution of CO ^ into the 
atmosphere (Rubey, 1951). 

Disposal of the hypothetical solid products is not difficult. The 

22 

graphite produced, 4.61 x 10 moles, would form a sphere about 390 km. in 

radius, or a layer on the surface of the earth about 0.5 km. thick. However, 

if it constituted 1 per cent of the upper part of the earth, it could be 

dispersed in a layer 50 km. thick. If the reduced phase, R, is considered 

9 3 

as iron, a volume of 1.03 x 10 km. is produced, corresponding to a sphere 
about 6 25 km. in radius or to 1 per cent of the earth to a depth of about 
200 km. If temperatures above the melting point of iron are reached, the 
metal could segregate toward the center. Other reduced phases which could 
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develop would have slightly larger volumes, but could still be distributed as 
small percentages of the material in the upper part of the earth. 

Despite the oversimplified assumptions of this model and the uncertainties 
in the calculations, it does appear that all the "excess volatile" water and 
00^ now present could have originated from reaction of original methane with 
accessory minerals at moderate depths within the earth. Such a model has 
two significant advantages. First, no rigorous assumptions are needed as to 
the original mineralogy and composition of the earth; the reactions involve 
only small amounts of accessory minerals, and gross mineral ogical characteristics 
are unaffected. Second, there is no necessity for equilibrium between the 
atmosphere and the surface rock at any time. 

The high calculated surface pressure for original methane, about 140 b., 
suggests that it did not form an atmosphere prior to reaction. The required 
methane could have been contained within the upper levels of the earth, 
although, if temperatures above 1000°C existed in this region soon after 
accretion (MacDonald, 1959), methane would not remain in equilibrium with 
magnetite-graphite assemblages, and conversion to water and 00 ^ would have 
proceeded rapidly. The product CO^ would also have exerted an unreasonably 
high surface pressure, suggesting that it was evolved gradually until formation 
of carbonate-bearing sediments could occur (Rubey, 1951). 

This mechanism is in general agreement with models of atmosphere develop- 
ment proposed earlier (Rubey, 1951; Holland, 1962). However, consideration 
of possible heterogeneous reactions at depth removes the necessity for an 
original methane-rich atmosphere. The model leads directly to an atmosphere 
enriched in water, CO and N 2 which resembles the second stage of Holland’s 
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(1962) model or the original atmosphere proposed by Rubey (1951). 

The materials-balance method used here gives no indication of the 

time required to produce the water and CO^ or of whether the rate of 

production has been constant throughout geologic time. Rubey (1951, p. 1138) 

has pointed out that the present rate of introduction of volatiles to the 

surface could product the required amounts during geologic time. There is 

thus no reason to require large production of volatiles early in the earth’s 

history, but the mechanism described would be consistent with either model . 

If temperatures in the atmosphere are too low for establishment of 

equilibrium, further compositional changes will depend largely on the 

initiation of photosynthesis to produce the present oxygenated atmosphere. 

The date for the origin of life and the magnitude of photosynthetic effects 

are uncertain; the discovery of apparently photosynthetic algae in Precambrian 

rocks (Rutten, 1962; Barghoorn and Tyler, 1965) indicates that photosynthesis 

9 

was operative, at least locally, as much as 2 x 10 years ago, a result whLch 

* : 

must be included in any consideration of the evolution of the atmosphere^ 
(Cloud, 1965). 


The system studied, including only C, H, and 0, gives no indication 

of the behavior of during oxidation of the methane. Higher values of 

P 0 promote the decomposition of NH^ to N 2 and , and it is probable that 

such reaction, occurring in parallel with the oxidation of methane, could 

produce the observed in the atmosphere. The amount of ’’excess volatile” 

21 20 

N 2 (Rubey, 1951, p. 1116) is 4.2 x 10 gm. or 1.5 x 10 moles, about 
0.002 of the amount of water and CO^ required. 
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GENERAL CONCLUSIONS 

(1) The study of equilibrium between the mineral graphite and a 

coexisting gas phase in the system C-H-G is of extreme significance in 

estimating conditions of formation of graphite-bearing mineral assemblages 

which occur in metamorphosed sediments, in certain igneous rocks, and in 

several classes of iron and stony meteorites. Ideally, the system graphite + 

gas is trivariant; in the present study, values of T, anc * are 

independently specified. Values of P , P P , P , and P in the 

GO 2 GO H. 2 O H. 2 C n.^ 

coexisting gas phase are then explicitly calculated with a high-speed computer. 
The numerical results support some general conclusions about the composition 
of the gas phase in equilibrium with graphite-bearing assemblages over a 
range of temperatures and gas pressures attainable in geological environments. 

(2) The contribution of methane, CH^, to the gas phase in equilibrium 

with graphite-bearing assemblages is significant over a wide range of 

geologically realizable conditions. In a moderately reduced gas phase, methane 

may constitute 30 to 50 per cent of the total gas pressure, with the excess 

composed chiefly of water and 00 ? . Under conditions of even lower p , the 

°2 

gas phase may be almost entirely methane. The existence of significant amounts 

of methane in the gas phase is indicated by the presence of minerals stable 

at low P , such as fayalite, wtistite, and native iron. 

°2 

(3) Graphite in a sediment will act to produce increasingly reduced 
mineral assemblages during rising temperatures associated with progressive 
metamorphism. Such a mechanism may explain the commonly observed reduction 
of original sediments with increasing metamorphic grade in both regional and 
contact-metamorphic environments. In particular, the existence of graphite 



or equivalent organic matter will reduce original hematite to magnetite at 
very low temperatures. 

(4) Because of the low gram-molecular weight of graphite, trace amounts 

will be sufficient to act as an effective internal oxygen buffer in systems 

as small as individual rock layers. Observed differences in P n in metamorphic 

U 2 

rocks over short distances may be conditioned by the presence or absence of 
graphite in the assemblage; a definite distinction between the two cases is 
necessary for a detailed study of such rocks. 

(5) The equilibrium relations suggest a model for the development and 
evolution of the earth's atmosphere in which original methane is oxidized 
to water and 00 ^ by reaction with trace amounts of graphite and ferrous- 
ferric minerals such as magnetite in the upper part of the primordial earth. 
Even for maximum required water and GO^, the materials needed are not 
excessive and could reasonably have been contained within a variety of possible 
primordial earths. Such a mechanism is independent of the gross composition 
and mineralogy of the earth and does not require high surface temperatures 

or establishment of equilibrium between atmosphere and surface. The reactions 
can therefore be applied to a wide range of proposed origins for the earth, 
its atmosphere, and its oceans. Possible reaction products are graphite, 
hydrogen, and iron, all of which could be incorporated into the primordial 
earth without difficulty. The time required for completion of the process 
is indeterminate, although the calculated high surface pressures of both 
reactant methane and product 00^ favor the interpretation that a relatively 
long period of time was involved (Rubey, 1951; Holland, 1962). However, 
the process could have been completed early in the Precambrian before most of 
the available sedimentary record originated. 
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In summary, graphite is a common and widespread accessory mineral in 
both terrestrial and extraterrestrial rocks. It is of particular importance 
because of its relationship to the composition of a coexisting gas phase in 
the system C-H-O. The presence of graphite in a mineral assemblage thus 
provides an additional constraint on the conditions of temperature, pressure, 
and composition in which the assemblage formed. Furthermore, changes in 
composition of the coexisting gas phase in response to changes in temperature 
have considerable significance for studies of progressive metamorphism and 
for speculations about the origin and development of the terrestrial 
atmosphere. For these reasons, the detection and identification of even 
trace amounts of graphite in a mineral assemblage are of particular importance 
to petrologic studies. 
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FIGURE CAPTIONS 

1. Calculated composition of a gas phase in equilibrium with graphite as 

a function of f for a temperature of 327°C (600°K) and total gas 
°2 

pressures of 1 bar, 1000 bars, and 10 kilobars. Values of Log.^ P^ for 

00 2 , 00, H 2 0, H 2 , and CH 4 are plotted against Log 1Q f Q ; dashed extensions 

of the curves indicate impossible values for which p co * p a> >! W 

Vertical lines above the curves indicate f values for solid-phase 

°2 

oxygen buffers. 

2. Calculated composition of a gas phase in equilibrium with graphite as a 

function of f for a temperature of 627°C (900°K) and total gas pressures 
U 2 

of 1 bar, 1000 bars, and 10 kilobars. Values of Log.^ P^ for C0 2> 00, 

H 2 0 , H 2 , and CH^ are plotted against Log^ f^ ; dashed extensions of the 

curves indicate impossible values for which P„ + P„ > P . Vertical 
* 00 2 00 gas 

lines above the curves indicate f values for solid-phase oxygen buffers. 

°2 ’ 

3. Calculated composition of a gas phase in equilibrium with graphite as a 

function of f for a temperature of 927°c (1200°K) and total gas pressures 
°2 

of 1 bar, 1000 bars, and 10 kilobars. Values of Log^ Q P^ for 00 2> 00, 

H 2 0 , ^2 anc * ^4 are Pl° ttec * against Log.^ f^ ; dashed extensions of the 

curves indicate impossible values for which P__ + P > P . Vertical 

v 00 2 CD gars 

lines above the curves indicate f values for solid-phase oxygen buffers. 

''i 
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4. Calculated composition of a gas phase in equilibrium with graphite 

as a function of temperature for f values specified by the 

U 2 

quartz-f ayal ite-magnetite buffer. Total gas pressure equals 1000 
bars (left) and 10 kilobars (right). Values of Log^ for CD 2> 

H 2 0 , ^2’ anc ^ are against temperature in the upper part of 

the diagrams; corresponding values of Log^ f^ are plotted in the 
lower part. 


5. Calculated composition of a gas phase in equilibrium with graphite as 
a function of temperature for f^ values specified by the combined 
magnetite- iron and magnetite-wllstite buffers. Total gas pressure 
equals 1000 bars. Values of L°g 10 P^ for 00 2 , CD, H 2 0 , H 2 , and CH^ 
are plotted against temperature in the upper part of the diagram; 
corresponding values of Log^ f^ are plotted in the lower part. 


Calculated composition of a gas phase in equilibrium with graphite as 

a function of temperature for fg^ values controlled by the combined 

magnetite-iron and wdstite-iron buffers. Total gas pressure equals 

1000 bars. Values of Log^ for OOj, CD, HjO, H 2 , and CH^ are plotted 

against temperature in the upper part of the diagram; corresponding 

f n values are plotted in the lower part. 

2 


7. Calculated mole percentages of methane, CH^, as a function of temperature 

in a gas phase coexisting with graphite at f values established by 

0 2 

various solid-phase oxygen buffers. Total gas pressure equals 1 bar. 
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8. Calculated mole percentages of methane, CH^, as a function of temperature 

in a gas phase coexisting with graphite at f values established by 

°2 

various solid-phase oxygen buffers. Total gas pressure equals 1000 bars. 


9. Calculated partial pressures of (X> 2 , CO, H^O , and CH^ as a function 
of temperature for gas phases coexisting with graphite. Total gas 
pressure equals 1 bar. The curves correspond to values of fg specified 
by the quart z-fayal ite-magnetite buffer (solid lines) and by the 
combined magnetite-iron and magnetite-wUst ite buffers (dashed lines). 


10. Calculated partial pressures of 00^, 00, H^O , , and CH^ as a function 

of temperature for gas phases coexisting with graphite. Total gas 
pressure equals 1000 bars. The curves correspond to values of f 

U 2 

specified by the quartz-fayal ite-magnetite buffer (solid lines) and 
by the combined magnetite-iron and magnet ite-wUst ite buffers (dashed 
lines) . 


11. Calculated temperature - total gas pressure curves which represent the 

upper limits of stability of various solid-phase buffer assemblages with 

graphite and a gas phase. The curves are actually projections of 

univariant curves along the f_ axis onto the P -T plane. 

0 2 gas 


12. Calculated values of P and P in a gas phase in equilibrium with 

H 2 0 00 2 

graphite at a total gas pressure of 100 bars. At each temperature, 
permissible compositions lie on curved isotherms (solid lines). Dashed 
lines indicate estimated compositions for a gas phase in equilibrium 
with graphite and various solid-phase buffer assemblages. 
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13. Calculated values of P and P in a gas phase in equilibrium with 

H 2 ° 00 2 

graphite at a total gas pressure of 1000 bars. At each temperature, 
permissible compositions lie on curved isotherms (solid lines). 

Dashed lines indicate estimated compositions for a gas phase in 
equilibrium with graphite and various solid-phase buffer assemblages. 

14. Values of f established in a gas phase in equilibrium with graphite 

U 2 

as a function of temperature for total gas pressures of 1 bar and 
2000 bars. Values for a system in which the H/0 ratio of the gas 
phase is fixed at 2/1 (heavy solid lines), corresponding to progressive 
metamorphism of a graphite-bearing sediment containing interstitial 
water, fall close to those calculated for the system C-0 (heavy dashed 
lines) (French and Eugster, 1965). Both sets of curves intersect 
progressively reduced solid-phase buffers (light solid lines) with 
increasing temperature. 

15. Calculated composition of a gas phase in equilibrium with graphite as 

a function of temperature for a system in which the H/0 ratio of the 

gas phase is maintained at 2/1. For temperatures below 1000°C, H 2 0 

dominates, and 00 2 and CH^ are present in approximately equal amounts. 

The f -T curve for this system (heavy solid line in lower part of 
°2 

figure) intersects more reduced buffers; with increasing temperature, 
an original magnetite-quartz-graphite rock should progressively 
develop fayalite, wUstite, and iron. 
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